clc
clear

data = importdata('..\run\check_halo.txt');

nx = size(data,2);
ny = nx;
nVar = 6;

var = zeros(nVar,nx,ny);

count = 0;
for iVar = 1:nVar
    for j = 1:ny
        count = count + 1;
        var(iVar,:,j) = data(count,:);
    end
end

figure
iVar = 3;
var_plt = squeeze(var(iVar,:,:))';
plt = pcolor(var_plt);
colormap(jet)
set(plt,'edgeColor','none')
% shading interp